Phase Separation in the Two-Dimensional Hubbard Model 
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By making use of known exact results and symmetry properties for the one-band Hubbard model, 
we show in a somewhat exact manner that there is no phase separation on a square lattice at 
arbitrary fillings at finite temperature for both attractive and repulsive on-site Coulomb interaction. 
This result is consistent with the quantum Monto Carlo calculation. 
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The considerable experimental evidence shows that 
La,2CuOi+s has a regime in which phase separation oc- 
curs 0. Such a separation appears between a phase 
very close to LaiCuO^ (i.e., an oxygen-poor phase) and 
an oxygen-rich phase that becomes superconducting at 
T c ~ 38-R" . Later on, this phenomenon was actively in- 
vestigated and confirmed in many other high Tc super- 
conductors (|]. People therefore believe that the pres- 
ence of phase separation is of essential importance for 
understanding the physics of the cuprate superconduc- 
tors [pi. Apart from the experimental works, there are 
also some theoretical investigations devoting to this sub- 
ject, see Refs. [^[Q for excellent reviews. As most of 
theoretical studies, being based on the planar t-J and 
Hubbard models, are numerical or approximate, the ob- 
tained results still remain controversial 0,0. Quite re- 
cently, Au et al. ||, using the symmetry properties and 
known exact results obtained by Lieb |(| , Kubo and Kishi 
JjJ for the Hubbard model, got an exact result on phase 
separation in the Hubbard model on bipartite lattices. 
They asserted that there is no phase separation at low 
dopings at any temperature for the repulsive Hubbard 
model. This exact consequence clarifies, to some extent, 
the existing controversy in the Hubbard model and in 
the t-J model with small values of J/t. However, sev- 
eral issues concerning this problem still remain to be ad- 
dressed in the sense that the exact results are sparse. For 
instance, if phase separation exactly occurs in the cases 
with moderate or even high dopings or in the attractive 
(negative-U) Hubbard model on symmetric bipartite lat- 
tices (e.g. on square lattices or cubic lattices), is still 
inconclusive. Since there have been numerical results in- 
dicating no phase separation for the one-band Hubbard 
model on square lattices at any fillings ||0], it would be 
quite interesting to seek more exact evidences supporting 
this observation. 

In this note, making use of known exact results ob- 
tained by Ghosh and symmetry properties for the 
one-band Hubbard model on a square lattice, we shall 
show that there is no phase separation at any fillings at 
finite temperature for both attractive and repulsive on- 
site Coulomb interaction. This statement, being exact, 
is consistent with the quantum Monto Carlo calculation 



We start from the one-band Hubbard model in an ex- 
ternal field ft on a square lattice with M sites. The Hamil- 
tonian reads 

H = -t ( c L c j° + h - c -) + U Y1 " l ~r™a 
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^(n iT + mi) - - J2(n^ - mi), 
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where the notations are standard. The foregoing discus- 
sion is independent of the sign of U, and is only valid 
for finite temperatures and for the finite chemical poten- 
tial. In the following, we shall first investigate the phase 
separation near half-filling, and then discuss the problem 
away from half-filling. The advantage of our method is 
the fact that we can treat the problem in a unifying way 
for both U > and U < 0. 

The method which we shall adopt for the case near 
half-filling, being based on the well-known particle-hole 
transformation, is very similar to that exploited in Ref. 
II, but we do not invoke any a priori assumption. For 
the Hamiltonian (1) defined on a square lattice, Ghosh 
M obtained, for small h, an exact upper bound 



\m(h,T)\ < 



const. 



T? I In I 



(2) 



where m(h,T) = jyEiv 5 ?}, the magnetization per site, 
and (• • •) is the thermal average over a grand canonical 
ensemble. As ft. 0, \m(h,T)\ — ► means the absence 
of spontaneous magnetic long-range order, implying that 
the system could exhibit paramagnetic behaviors at tem- 
perature T > 0. For small ft, inequality (2) ensures 
the analyticity (only in the sense that the first deriva- 
tive exists) of m(ft,T). Therefore, we can conclude that 
m(ft, T > 0) is continuous and analytic in the neighbor- 
hood of ft = 0. 

We now apply the unitary particle-hole transforma- 
tion, 



(3) 



with e(i) = —1 for i G the one sublattice and —1 for i G 
the other sublattice, to the Hamiltonian (1), and obtain 



1 



H' = -t (cl c ^+ h - c -) + u 'J2 n ^ ni i 



ti 



with 



-V'^2(ni1 +nn) - -^^{nti - mi), 



, h-U 
U' = -U, , h' = 2fi-U, 



(4) 



(5) 



where a constant term is dropped. Eq. (2) is transformed 
into 



< 1 



const , . , 

TV2 |ln|C/'-2//||V2 - 

const. 1 



1 



T 1 / 2 |ln|t/'-2/i'|| 1 /2 



(6) 



the 



t J 



(or pseudospin) operators PJlT|] as 77+ = J2i e (*) c 
f = (^ + ) t ,and?7 z = t}(N—M) with JV = £i(n iT +n a j. 
They obey the usual SU{2) Lie algebra. We from now 



on set h = in Eq.(l). As a matter of fact, , H] = 
2(p c — p)n~ , where p c — ^. Using this commutator and 
the cyclicity under the trace we obtain [TL2[ 



F(p, T){1 - exp[2/3( Mc -/*)]} = T) - 1, 



(7) 



for small |J7' — 2/i'|, where p = jjJ2i( n 
electron density per site. Inequality (6) is of basic im- 
portance. For the system described by the one-band 
Hubbard model (4) on the square lattice, the density 
p(p' , T > 0) is thus the continuous function of the chem- 
ical potential p' near i.e., at small doping (note that 
at half-filling, p! = ^-). This can be justified by noting pf^ T) = 1 
the fact lim^o p(p' = U'/2 + 6, T > 0) = lim^o p(p' = 
U' /2 — S, T > 0). An alternative criterion for phase sepa- 
ration is based on this fact Q : if a discontinuity is found 
in p(p',T) as a function of p' , then the densities inside 
the gap are unstable, giving rise to a phase-separated 
state, and if it is not found, then no phase separation oc- 
curs. Actually, phase separation falls in the class of the 
first-order phase transition, as the first-order phase tran- 
sition into two phases with different densities is featured 
by the discontinuity of p(//,T). Thus the continuity of 
p(p',T > 0) near U'/2 suggests that the one-band Hub- 
bard model on the square lattice at T > does not ex- 
hibit phase separation at small doping (near half- filling) , 
which is consistent with the quantum Monto Carlo cal- 
culation IIU] and with the result of high-temperature 
expansion for the two dimensional t-J model with small 
values of J/t [[[o| which can be viewed as the strongly 
coupling limit of the Hubbard model. As in the begin- 
ning we do not fix the sign of U, this exact result re- 
mains true for both negative and positive U' as well as 
the vanishing external field [hi = 0). If the system is 
indeed paramagnetic, as being plausibly expected, then 
this result could be extended to moderate dopings. One 
may observe that this approach only works for the two 
dimensions (or the one dimension) but not for the three 
dimensions, as the Ghosh's results were obtained only for 
the low dimensional cases. 

Now let us look at this problem when the system is 
doped away from half-filling. We introduce the 77 pairing 



with F(p,T) = jj{n + n~) , and (3 the inverse temper- 
ature. Since the right-hand side of (7) is an intensive 
quantity, F(p,T), being obviously non-negative, should 
also be intensive. Furthermore, it is well-known jl3| that, 
based on the particle-hole transformation, p(p, T) = 1 
(half-filling) as fi = p, c and vice versa. Consequently, 
combining this result and Eq.(7) we get the following re- 
lation 

p(p, T) > 1 as p > pa p(p, T) = 1 as p = p c , 
and p{p,T) < 1 as p < p c . (8) 

This exact constraint may have some implications. First, 
a direct consequence is F(p,T) > for p 7^ p c . Second, 
since p(p, T) is analytic in the neighborhood of p = p c 
(ensured by (6)), at small doping, we can expand p(p, T) 
in powers of (p — p c ), obtaining 



dp 



\^{P - Pc) + 0[(p - p c )\ (9) 



where we have used p{p — p c -T) — 1. Up to the second 
order in (p — p c ), we see that ' 9p g^' T " > U=Mc = finite > 
(otherwise it contradicts (6) and (8)). 

Since the functions p(p, T) and F(p, T) are closely re- 
lated, let us now study the properties of the latter. From 
Eq.(7) we know that F(p,T) is finite, and comply 



0<F(p,T) < 



1 



l-exp[2/3(Mc-M)]| 



(10) 



for p ^ p c . Furthermore, in accordance with the defini- 
tion of F{p, T), we have 



F(P,T) 



1 Tr[e~ l3( - Ho ~ tlNS, n + ri~ 



Tre -f3(H -vN) 



(11) 



where we have denoted Eq.(l) by H = Hq — pN (Recall 
that h = 0). Employing the following symmetric particle- 
hole transformation, 



"it 



e(i)c 



e(*) c iij 



(12) 



we find H -> H = H - (U - p)N + (U - 2p)M, and 
r/ + r/~ — > rj~r} + . By applying the unitary transformation, 
Eq.(12), to F(p,T), we obtain 



F(p,T) = 



1 Tr{e 



-/3[H -(U-n)N] 



n n + } 



M 



= F(2p c -p,T} 



Tre -l3{H ~(U-n)N] 
,2P(ji-/i c ) 



(13) 



where we have used n~ri + = n + r]~ — 2n z , the definition 
of F(p,T), as well as Eq.(7). Eq.(13) completely deter- 
mines the form of F(p,T). As p! = p c — p, it becomes 



2 



F(p c - p,T) = F(p c + p,T)exp(-2f3p), which reflects 
the symmetry of F(p, T) as a function of p. Similarly, 
we can obtain, after operating the transformation (12) to 
p(p,T), the following equation 



p(p,T) = 2-p(2p c -p,T). 



(14) 



One may verify that Eqs.(7), (13) and (14) are self- 
consistent. 

By differentiating Eq.(ll) with respect to p, and us- 
ing the unitary particle-hole transformation (12) to the 
thermal averages involved, and then noting Eqs.(7) and 
(13), one can prove exactly the following expression 



dF(p,T) dF(2p c -p,T) 



dp 



-^«iV 2 ) - {N) 2 h^ T , 



(15) 



where (• • -^t means the thermal average with the chem- 
ical potential p at temperature T. Considering (iV 2 ) — 
(iV) 2 > 0, we get > d -E^^H . By (14) and 

(10), we obtain the inequality 



dF(ji,T) < 2(ie 2f3{ -^-^ 



dp. 



( e 20(^c-/i) - 1)5 



(16) 



for p < p c . 

Till now we have already had some basic knowledge 
about the function F(p,T), i.e., it should satisfy (7), 
(10), (13), (15) and (16) simultaneously. Under these 
conditions, we can solve equation (13) exactly, and the 
solution, being surprisingly simple, is given by 



F(p,T) = 



C(T) 



1 + e 2 / 3 (/ J c-M) ' 



(17) 



where C(T) is a positive, finite constant at temperature 
T, satisfying 



< C{T) < 1. 



(18) 



We like to mention here that we can not obtain the closed 
form of C(T) for the time being, but such a form is 
enough for our purpose. Consequently, by Eq.(7), we 
have 



p{p, T) = 1 + C(T) tanh[/3(/x - p c )}. 



(19) 



One may check that the solutions given by Eqs.(17) and 
(19) satisfy all aforementioned properties of functions 
F(/i,T) and p{p,T). 

With these facts, we are ready to address the problem 
of phase separation away from half-filling. As we have 
obtained an exact, explicit expression for p(p 7 T), say, 
Eq.(19), it is obvious that p(p,T) as a function of p is 
continuous at any fillings. This result shows that there is 
no occurrence of the first-order phase transition, which in 



turn leads to the conclusion that the one-band Hubbard 
model can not exhibit the phenomenon of phase separa- 
tion on a square lattice at finite temperature, consistent 
with the numerical calculations . 

A few remarks concerning Eqs.(17) and (19) are in or- 
der. (1) As one may note that, Eq.(19) looks to be rea- 
sonable, as it comes directly from the symmetry of the 
Hubbard model, and is well in agreement with the in- 
equality (6). Moreover, p(p,T), given by (19), has the 
qualitatively similar behaviors compared with numeri- 
cal results on square lattices with small sizes, though a 
quantitative comparison is impossible due to uncertainty 
of C(T) as well as the lack of accurate data. Clearly, 
Eq.(19) remains valid in the thermodynamic limit. We 
like to mention here that we can not prove the uniqueness 
of the solution Eq.(17), but the other solutions, if exist, 
might have the similar forms as Eq.(17) [jl4|, which would 
not affect our conclusion. (2) As indicated in Eq.(19), 
p(p, T) might have different behaviors for positive and 
negative U because of possible p c — dependence of C(T). 
(3) We emphasize once again that Eqs.(17) and (19) work 
only for finite temperatures and finite chemical poten- 
tials. One can not extract any useful information for 
the zero-temperature case from this study. Based upon 
the present result, however, we could say that it may 
be inappropriate to use the one-band Hubbard model as 
a model to explain the phenomena of phase separation 
observed in high Tc superconductors, at least at finite 
temperatures. 

In summary, we show in a somewhat exact way that, 
using the known exact results and symmetry properties 
for the one-band Hubbard model, there is no phase sep- 
aration at any fillings on a square lattice at finite tem- 
perature for both negative and positive U. This result 
is consistent with the quantum Monto Carlo calculation 
and some analytic result, and is also compatible with the 
recent exact result obtained in Ref. Nevertheless, we 
like to mention that our result does not cover the case 
at zero temperature. As there are also strong numerical 
evidences showing no occurrence of phase separation at 
zero temperature, how to obtain an exact proof for both 
U > and U < is still a fascinating topic. 
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